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ABSTRACT 

Context. Establishing the number of faint AGNs at 2 ; = 4 — 6 is crucial to understand their cosmological importance as 
main contributors to the reionization of the Universe. 

Aims. In order to derive the AGN contribution to the cosmological ionizing emissivity we have selected faint AGN 
candidates at 2 > 4 in the CANDELS GOODS-South field which is one of the deepest fields with extensive multiwave¬ 
length coverage from Ghandra, HST, Spitzer and various groundbased telescopes. 

Methods. We have adopted a relatively novel criterion. As a first step high redshift galaxies are selected in the NIR H 
band down to very faint levels {H < 27) using reliable photometric redshifts. This corresponds at 2 > 4 to a selection 
criterion based on the galaxy rest-frame UV flux. AGN candidates are then picked up from this parent sample if they 
show X-ray fluxes above a threshold of Fx ~ 1.5 x 10“^^ erg cm“^ s”*^ (0.5-2 keV), correponding to a probability of 
spurious detections of 2 x 10”"* in the deep X-ray 4 Msec Chandra image. 

Results. We have found 22 AGN candidates at 2 > 4 and we have derived the first estimate of the UV luminosity 
function in the redshift interval 4 < 2 < 6.5 and absolute magnitude interval —22.5 < M 1450 < —18.5 typical of local 
Seyfert galaxies. The faint end of the derived luminosity function is about two/four magnitudes fainter at 2 ~ 4 —6 than 
that derived from previous UV surveys. We have then estimated ionizing emissivities and hydrogen photoionization 
rates in the same redshift interval under reasonable assumptions and after discussion of possible caveats, the most 
important being the large uncertainties involved in the estimate of photometric redshift for sources with featureless, 
almost power-law SEDs and/or low average escape fraction of ionizing photons from the AGN host galaxies. Indeed 
both effects could in principle significantly reduce the estimated average volume densities and/or ionizing emissivities 
especially at the highest redshifts. 

Conclusions. We argue that, under reasonable evaluations of possible biases, the probed AGN population can produce at 
2 = 4 — 6.5 photoionization rates consistent with that required to keep highly ionized the intergalactic medium observed 
in the Lyman-a forest of high redshift QSO spectra, providing an important contribution to the cosmic reionization. 

Key words, galaxies - AGNs - reionization 


1. Introduction 

The process of cosmic reionization remains one of the most 
important and puzzling problems for cosmology. The reion¬ 
ization of the universe, after its neutral phase followed by 
recombination, appears well completed at 2 ~ 6 as in¬ 
ferred from the high ionization level of the intergalactic 
medium (IGM) obser ved in the Lyman -a forest of high 
redshift QSO spectra (Fan et al. (2006) I. The epoch when 
the process becomes important is still unknown although 
some hints come from the level of thompson scattering of 
the CMB radiation observed by WMAP and Planck and 
by the analysis of the kinematic Sunayev-Zeldovich effect 


(Hinshaw et al. (2013) George et al. (2014)). The common 
wisdom is that in the epoch corresponding to the redshift 
interval 2 ^ 6 — 10 the Universe has quickly changed its 
thermal and ionization state. However the identification of 
the source populations emitting enough UV radiation to 
reionize the Universe is still an open issue. 

In more than fifteen years extensive studies have been 
performed to identify the ionizing population among the 
high redshift star forming galaxies or the bright AGN 
population. Since the apparent number density of bright 
QSOs and AGNs is rapidly decreasing at 2 > 3 (e.g. 
Masters et al. (2012) I it is usually assumed that QSOs do 
not produce enough ionizing emissivity at 2 > 4. Indeed 
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estimates based on simple parametric extrapolations both 
in luminosity and redshifts of the AGN luminosity func¬ 


tion imply values for the cosmological emissivity 10^^ ergs 


Mpc ^ Hz ^ s ^ (e.g. Haardt & Madau (2012) |. This is 


one order of magnitude smaller than what required to 
keep the intergalactic medium highly ionized at z ^ 6 de¬ 
spite the fact that each bright QSO is able to ionize its 
neighbourhood u p to a distance of the order of e.g. 10 
Mpc at -z ^ 5.5 ( Prochaska, Worseck, &: O ’Meara (2009)[ 


Songaila fc Cowie (2010) ' Worseck et al. (2014)p. 

As a consequence, the more common star forming galax¬ 
ies are thought to be responsible for the cosmic reioniza¬ 
tion. Their contribution depends on their abundance at 
low luminosities and on the average fraction of ionizing 
UV flux escaping from each galaxy into its neighbourhood. 
At very high redshifts, z > 5, direct measures of the ion¬ 
izing Lyman continuum flux from galaxies are not viable 
due to the saturated hydrogen absorption by the IGM ob¬ 
served in QSO spectra. Estimates at lower z ^ 3 —4 become 
feasible but the values so far appear discrepant relying on 
specific assumptions and interpretations. For example, re¬ 
cent analyses suggest average ion i zing escape fractions o f 
10 — 15% (Nestor et al. (2011) Mostardi et al. (2013) I. 
But their z ^ 3 candidates show ionizing continuum of¬ 
ten slightly shifted in position respect to the non-ionizing 
flux, casting serious doubt of significant contamination by 
foreground low redshift interlopers (Vanzella et al. (2012) |. 
Other estimates based on spectroscopic and very deep 
broad or narrow band imaging gave only uppe r limits < 5% 
( Vanzella et al. (2010)[ Boutsia et al. (2011)). At present 
the results on the estimated ionizing escape fraction for 
this relatively bright star forming population at z ~ 3 are 
not conclusive. 

Faintest dwarf high redshift galaxies could behave in 
a different way showi ng larger escape fractions (see e.g. 
Fontanot et al. (2014)) and could be significant contrib¬ 
utors to the cosm i c reionization ([Ferrara fc Loeb (2013)[ 
Wise et al. (2014')| |Yue et al. (20l4)| . However some as- 
sumptions on the abundance at z '^ 7 — 10 of the faintest 
galaxies down to Mj/y ^ —13 and on their escape frac¬ 
tions (^ 15 — 80%) are needed to make the star forming 
gala xy population a major contributor to the reioniza tion 
(e.g. Finkelstein et al. (2012) Robertson et al. (2013)). A 


at z > 3 (|Glikman et al. (2011) Civano et al. (2011), 
Fiore et al. (2012)[ ) showing the presence of a consider¬ 
able number of previously unknown faint AGNs able to 
produce a rather steep luminosity function. The pres¬ 
ence of this faint population is changing our estimate of 
the AGN contribution to the ionizing UV background al¬ 
though the selection of faint AGNs at the highest red¬ 
shifts z ~ 6 is difficult with the current instrumenta¬ 
tion. Most important, the inclusion of X-ray detection in 
the selection methods enables to extend the knowledge 
of the luminosity function to even fain ter X-ray limits 
(Givano et al. (2011) Fiore et al. (2012)1 of the order of 
L(2-10 KeV) > 10'^'^ erg s“^ at z > 4. In particular, the 
novel technique proposed by Fiore et al. (2012) relies on 
the selection of faint high redshift AGNs among the high z 
galaxies selected in NIR images (e.g. H band) which show 
any detection in deep and high resolution, Ghandra X-ray 
images. 

A first estimate of the volume density of faint 
AGNs could also have important implications concern¬ 
ing the abu ndance and mass of the supermassive black 
hole seeds (Volonteri (2012) ) and their early growth 
(Volonteri & Silk (2014)p. Important constraints are also 
expected for the initial spectrum of the density perturba¬ 
tions rel ated to the nature, co ld or warm, of the dark matter 
(see e.g. Mend et al. (2013)D. 

In |Giallongo et al. (2012) we made a first attempt to 
estimate the possible contribution to the cosmic reioniza¬ 
tion by faint high z AGNs based on the results provided by 
deep X-ray surveys and on a joint galaxy-AGN model pre¬ 
diction. We found that faint AGNs could be able to reionize 
the high redshift IGM provided that their average volume 
density only gradually decline from z^3toz^6.5as 
predicted by the model. 

In the present paper we have produced a sam¬ 
ple of galaxies selected in the Gandels/GOODS- 
S/GhandraDeepField-South field using one among the 
deepest NIR, optical and X-ray database and we have built 
a sample of AGN candidates based on X-ray detection in 
the same H band galaxy position, following the recipe as in 


recent evaluation by I'inkelstein et al. (2014) | indicates es¬ 
cape fractions ^ 13% to keep the IGM ionized up to z ~ 6 
and midly ionized to z ^ 8 although uncertainties are very 
large. Specific searches are starting in lensing galaxy clus¬ 
ters to measure the escape of Lyman continuum flux in 
intrinsically faint galaxies which are gravitationally lensed 
and amplified in luminosity by the intervening clu ster (see 
e.g. Vanzella et al. (2012b) Ishigaki et al. (2014)1. 


All these uncertainties in the galaxy contribution leave 
open the possibility that AGNs could still play a leading 
role in contributing to the cosmological ionizing background 
although the estimate of their contribution is mainly af¬ 
fected by the poor knowledge of the faint end slope of 
their luminosity function especially at high redshifts due 
the lack of deep AGN surveys at various wavelengths 


Fiore et al. (2012)| In Sect. 2 we describe the GANDELS 
photometric parent catalog of high redshift galaxies. In 
Sect 3. we describe the AGN selection by means of the 
Ghandra 4Ms image. In Sect. 4 we present an estimate of 
the faint end AGN UV luminosity function. In Sect. 5 we 
show predictions about the AGN ionizing emissivity and 
photoionization rate and finally in Sect. 6 we discuss some 
possible caveats. Throughout the paper we adopt round 
cosmological parameters Da = 0.7, Dq = 0.3, and Hubble 
constant h = 0.7 in units of 100 km/s/Mpc. Apparent 
magnitudes are in the AB photometric system. 

2. The GANDELS GOODS-S catalog 


The _ GANDELS area (Groginetal. (2011) 


Koekemoer et al. (2011)) with deep Ghandra coverage 


(Shankar Sz Mathur (2007)). 

The traditional view based on an evolutionary 
ionizing QSO background increasing from the local 
value out to z = 2 — 3 and then quickly de¬ 
creasing (see e.g Haardt & Madau (2012)) is changing 

deep surveys 


Haardt & Madau (2012)) 
thanks to the recent multiwavelength 


is providing the optimal dataset in terms of deep NIR 
photometry and relatively wide area. Indeed the NIR 
images obtained with the Wide Field Camera 3 (WFC3) 
are particular important in terms of angular resolution 
and depth, sampling the AGN UV emission (< 3000 A) 
at z > 4. Among the various GANDELS fields we have 
selected in this preliminary analysis the GOODS-South 
field where the deepest 4Ms Chandra images are optimal 
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for the selection of high redshift AGNs. The field covers 
an area of ^170 sq. arcmin at a mean depth of H=27.5. 
The dataset also includes very deep imaging by the 
IRAC instrument on board the Spitzer Space Telescope 
coming from the S pitzer Extended Deep Survey (SEDS; 
Ashby et al. (2013)|) and covering the CANDELS fields 
to 26 AE mag (3 ct) at both 3.6 and 4.5 pm. The Spitzer 
infrared bands sample the rest-frame optical region of 
AGNs at z > 4. We adopted the official photometric 
catalog released by IGuo et al. 20131 selecting galaxies in 
the WFC3 H band. The total number of sources detected 
in the GOODS-South fields at a mean depth of H=27.5 is 
about 28500. 

Due to the complexity of the GOODS-South exposure 
map in the H band, the field is composed by several sub- 
areas with the associated magnitude limits. Simulations 
have been perfor med to assess the det ection complete¬ 
ness as described in Grazian et al. (2014) In the shallowest 


area, the incompleteness corrections as a function of mag¬ 
nitude have been estimated comparing the observed galaxy 
number counts with that in the deepest HUDF area. For 
the other areas of intermediate depths we scaled the correc¬ 
tion factors to fainter magnitudes taking into account the 
increasing s/n ratio. Table 1 provides the magnitude corre¬ 
sponding to a 50% level of completeness in each of the areas 
covered by different sub-regions. The numbers of galaxies 
at H <27 are also shown in each region and at different 
redshifts. 

Several catalogs of photometric redshifts have 
been provided for the CANDELS GOODS-S field. 
Dahlen el al. (2013) have made a comparison analysis 
among the various codes. Different galaxy templates have 
been used so that the codes were sufficiently different to 
provide independent redshift estimates. They have finally 
combined independent estimates obtained by using nine 
different codes. Several spectroscopic redshifts available 
in the observed field were compared to the resulting 
photometric redshifts of the same objects to derive the 
optimal solution for the overall sample. The resulting 
average uncertainty is | Az|/(1 -|- zspec) 0.03 at relatively 
bright magnitudes (iLieo < 24), which increases to 0.06 
at Hieo ^ 26 (Dahlen el al. (2013) I. The fraction of 
outliers (defined as differences |Az|/(l -I- zspec) > 0.15) 
is confined to ~ 3% and in particular is low at z > 3.7 
due to the Lyman continuum break signal. A Bayesian 
method has been developed to provide redshift probability 
distributions (PDFs) that can be useful for error estimates. 
The PDFs of the AGN candidates of the present sample 
are shown in the Appendix. 

In total, there are 1113 galaxies in the GOODS-South 
field with a robust spectroscopic or photometric redshift 
at z > 4 which are analysed for the search of any AGN 
activity. It is to note that the photometric redshifts have 
been estimated using theoretical spectral energy distribu¬ 
tions derived from stellar populations only. No contribu¬ 
tion from the non stellar light typical of the AGN popu¬ 
lation has been included. Nevertheless, since the selection 
of faint AGNs is confined to z > 4, the redshift estimate 
is based on the dropout of the SED due to the neutral HI 
absorption by the IGM intervening along the line of sight, 
a feature independent of the main intrinsic spectral proper¬ 
ties of the sources like e.g. continuum spectral shape or es¬ 
cape fraction of ionizing photons from the AGN host galaxy. 
Indeed a large UV dropout is expected in any case at z > 4 


by optically thick intervening absorption by the IGM even 
if the ionizing esca pe fraction from th e AGN host galaxy 
were high. Recently Hsu et al. (2014) have provided a cat¬ 
alog of photometric redshifts in the CANDELS GOODS-S 
area optimized for relatively bright X-ray detected AGNs. 
They add to the CANDELS photometric catalog new opti¬ 
cal intermediate-band filters and explore a mixing of galaxy 
and AGN templates. Their recipe appears to provide good 
results in terms of redshift accuracy and low fraction of out¬ 
liers when comparing with the spectroscopic redshift sub¬ 
sample. This is in particular true for bright sources with 
iL < 23. However for our much fainter sample having typ¬ 
ically iJ = 24 — 26 the use of shallower intermediate band 
filters and mixed AGN/galaxy templates does not improve 
the accuracy of photometric redshifts and the fraction of 
outliers. 


3. The AGN selection 

The selection of AGN candidates relies on the selection 
of galaxies at z > 4 having H <27 from the general 
CANDELS catalog as described in the previous section. 
This corresponds to select sources on the basis of their 
detected rest-frame UV luminosity at A 2000 — 3000 
A. Following the procedure proposed in Fiore et al. (2012) 
we searched for any sign of AGN activity just measur¬ 
ing the X-ray flux in the expected H band position of 
the CANDELS sources having photometric or spectroscopic 
redshifts z > 4. In the following we briefly summarize the 
main steps. Since X-ray data are structured in event files, 
preserving time, spectral and spatial information, a multi- 
di mensional source det ection technique has been developed 
by Fiore et al. (2012) using spatial, spectral, and timing 
information. In particular, clustering of X-ray events in 
energy, time and spatial coordinates was investigated to 
efficiently detect X-ray bursty sources and X-ray sources 
characterized by specific features like strong lines or sharp 
edges, by reducing the elapsed time and/or the energy in¬ 
terval (i.e. the background) where the source detection is 
performed. Spectra were added at each galaxy position in 
different apertures from three to nine pixels. Sources are se¬ 
lected having more than 12 background subtracted counts 
in the wide 0.3-4 keV band. The background is evaluated 
from a map obtained from the original image after the ex¬ 
clusion of all bright sources. The background for each source 
was computed by normalizing an average background at the 
source off-axis to the source spectrum integrated in the 7- 
11 keV energy band, where the source contribution is neg¬ 
ligible due to the decrease in the Ghandra effective area. 
Source extraction aperture and contiguous energy band 
have been selected to minimize the Poisson probability for 
background fluctuation and to maximize the signal-to-noise 
ratio. Extensive simulated data have been performed us¬ 
ing about 10® background spectra at the positions of the 
detected sources using the same exposure time, PSF etc. 
studying the number of spurious detections as a function of 
threshold, aperture and energy band width. A combination 
of parameters has been adopted to keep the number of spu¬ 
rious detections smaller than one over 5000 spectra corre¬ 
sponding to a probability of spurious detection of 2 x 10“"^. 
This implies that for the 1113 galaxies at z > 4 we ex¬ 
pect only 0.2 spurious AGN candidates in the field. Finally 
for each candidate position in the Chandra image we have 
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Table 1. Area and magnitude limits of the CANDELS GOODS-South and HUDF fields 


Field 

Area 

arcmin^ 

Hieo Mag. compl. 
50% 

< VI 

in 

■\rh,zgh, — z 
gal 

4 < z < 6.5 

GOODS-South #1 

11.05 

25.75 

801 

24 

GOODS-South #2 

25.03 

26.00 

2794 

77 

GOODS-South #3 

50.47 

26.50 

7984 

240 

GOODS-South #4 

77.18 

27.00 

15310 

698 

GOODS-South #5 (HUDF) 

5.18 

27.75 

1672 

74 

TOTAL 

168.91 

- 

26963 

1113 


Table 2. Photometric and spectroscopic properties of the 22 AGN candidates 


ID 

RA 

Dec 

zphot 

zspec 

G 

H 

Afl450 

log Ex 
erg/cm®/s 

log Lx 
erg/s 

A 

Previous Catalogs 

273 

53.1220463 

-27.9387409 

4.49 

4.762^ 

c 

23.96 

-21.37 

-15.97 

43.80 

#2 

M208,X403 

4285 

53.1664941 

-27.8716803 

4.28 


cf 

25.57 

-20.22 

-16.46 

42.90 

#3 

- 

4356 

53.1465968 

-27.8709872 

4.70 


cf 

26.36 

-18.44 

-16.38 

43.40 

#4 

M70437,L306,X485 

4952 

53.1605007 

-27.8649890 

4.32 


c 

25.47 

-20.20 

-16.50 

42.90 

#3 

- 

5375 

53.1026292 

-27.8606307 

4.41 


c 

25.16 

-20.16 

-16.66 

42.75 

#4 

X331 

5501 

53.1280240 

-27.8593930 

5.39 


c 

25.71 

-20.23 

-16.45 

43.10 

#4 

- 

8687 

53.0868634 

-27.8295859 

4.23 


c 

26.90 

-19.19 

-16.43 

42.90 

#4 

- 

8884 

53.1970699 

-27.8278566 

4.52 


c 

25.74 

-19.04 

-16.77 

42.65 

#4 

- 

9713 

53.1715890 

-27.8208052 

5.86 

5.70^ 

c 

26.54 

-19.87 

-16.46 

43.15 

#4 

HUDF322 

9945 

53.1619508 

-27.8190897 

4.34 

4.497® 

cd 

24.99 

-20.93 

-16.65 

42.75 

#4 

- 

11287 

53.0689924 

-27.8071692 

4.94 


c 

25.06 

-20.48 

-16.42 

43.10 

#4 

M8728 

12130 

53.1514304 

-27.7997601 

4.43 

4.62^ 

c 

25.54 

-20.60 

-16.58 

42.85 

#5 

HUDF3094 

14800 

53.0211735 

-27.7823645 

4.92 

4.823® 

c 

23.43 

-22.51 

-16.38 

43.10 

#3 

M10548 

16822 

53.1115637 

-27.7677714 

4.52 


c 

25.67 

-18.97 

-15.91 

43.85 

#4 

M70168,L245,X371 

19713 

53.1198898 

-27.7430349 

4.84 


c 

25.31 

-18.14 

-16.48 

43.00 

#4 

E1516,X392 

20765 

53.1583449 

-27.7334854 

5.23 


f 

24.44 

-21.06 

-16.29 

43.25 

#3 

E2551 

23757 

53.2036444 

-27.7143907 

4.13 


c 

24.56 

-20.72 

-16.49 

42.85 

# 1 

- 

28476 

53.0646867 

-27.8625539 

6.26 


f 

26.77 

-19.03 

-16.60 

43.10 

#4 

M70407 

29323 

53.0409764 

-27.8376619 

9.73 


cf 

26.33 

-19.50 

-15.96 

44.00 

#3 

M70340,L103,X156 

31334 

53.2131871 

-27.7816486 

4.73 


f 

26.41 

-19.60 

-15.69 

43.75 

#4 

- 

33073 

53.0547529 

-27.7368325 

4.98 


c 

26.89 

-19.19 

-16.44 

43.10 

#2 

E2199 

33160 

53.0062504 

-27.7340678 

6.06 


cf 

25.90 

-19.62 

-16.26 

43.65 

#3 

E2498,L57,X85 


ID from IGuo et al. 20131 catalog 
Photometric redshifts (zphot) from Dahlen el al. 


Rhoads et al. (2009) ' Balestra et al. (2010) 


catalog 


Spectroscopic redshift s after: qVanzeiia et al. (2008) | ^Hathi et al. (2008) [ '^M. Dickinson & D. Stern private communication, 


Column G indicates the compactness of the sources. c=compact, d=diffuse, f= faint object 

X-ray fluxes are co mputed in the 0.5-2 keV b and an d X-ray luminositie s in the 2-10 keV band M,E in Previous Gatalogs are from 
Fiore et al. (201^ L from Luo et al. (2008)] X from Xue et al. (2011) Column A indicates the GOODS-S field area according to 


table 1 


checked that the accuracy of the relative astrometry is < 1 
arcsec (see the Appendix for an example). 

We have detected 22 reliable AGN candidates with 
z > 4 whose positions, redshifts, X-ray (0.5-2 keV) fluxes 
and H magnitudes are shown in Table 2. Eleven of the de¬ 
tected sources were already in Eiore et al. (2012) catalog 
although with somewhat different photometric redshifts, 
four in the GOOD S-S catalogs of X -ray selected sourc es by 
Luo et al. (2008) and seven are in |Xue et al. (2011) (four 


in common with Luo et al. (2008 j[ ). h'ive redshifts have 
been confirmed spectroscopically either by the presence 
of strong Lyman-a intervening absorption trough (for the 
9713 and 12130 objects observed in slitless low resolution 
spectroscopy with ACS camera at HST) or by the detec¬ 
tion of Lyman-a emission. The detection of CIV at z > 4 
is usually not confirmed because out of the observed spec¬ 
tral region or expected in regions dominated by strong sky 
emission lines. The morphological appearance is indicated 
in column C where most candidates after visual inspection 


appear compact in H band with FWHM~ 0.2 — 0.4 arcsec. 
Only one (9945) has a complex morphology. The remaining 
are too faint for any classification. It is to note the larger 
number of AGN candidates emerging from the combined 
H band selection plus X-ray detection respect to the direct 
X-ray selection. X-ray luminosities have been computed for 
the brightest sources adopting a standard single power-law 
with fixed slope 1.8 of the power-law photon index and 
fitting a gas absorption. We derived hydrogen column den¬ 
sities of logNn - 23.5 for 273, 4356, 16822, 29323, 33160 
and logNff ^ 22.5 for 33073, 20765. For the remaining 
X-ray fainter objects we adopted logN^ ~ 22. The reli¬ 
ability of the estimated X-ray absorption derived for the 
estimate of X-ray fluxes and luminosities in z > 4 faint 
sources is affected by the higher rest-frame energy sam¬ 
pled, by the adopted slope in the power-law fitting and by 
other physical assumptions like covering factor and metal- 
licity of the absorbing gas, as discussed in the last section. 
The derived X-ray luminosities in table 2 have an aver- 
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Fig. 1. Image thumbs of two AGN candidates. From left to right the filters are Xo. 5 _ 2 , Xo. 8 - 4 , U, B 435 ,V 606 , I 775 , 2350 , 
Yi. 05 , Ji. 25 , Hi, 6 , IRAC 3 . 6 , 4 . 5 , 5 . 8 , 8- Left box: object 16822 at z = 4.5. Right box: object 28476 at z = 6.26. The size 
of the optical/IR thumbs is 9 x 6 arcsec^. Red circles in the X-ray images show the expected position from the H band 
detection with a circular size of 1 arcsec 


age value of log Lx ~ 43.2 erg s“^ and are more typical 
of Seyfert-like sources rather than of starburst galaxies al¬ 
though we can not exclude a significant contribution from 
stars to the X-ray lum inosity of some pec uliar source in 
our sample. Analysis by Fiore et al. (2012) of the brightest 
sources at 7 ; ~ 3 suggests that the possible presence of sig¬ 
nificant X-ray absorption is uncorrelated with the presence 
of Lyman-a and CIV emission lines and low dust extinc¬ 
tions as in the case where absorption is in general close 
to the emitting ionized region and inside the sublimation 
radius. 


Figure 1 shows as an example the thumbs of 2 sources 
with decreasing wavelength from the IR Spitzer images to 
the HST images as well as the X-ray image. The first source 
on the left, 16822, has a good signal-to-noise ratio and ap¬ 
pears compact in the H band. The second AGN candidate 
on the right, 28476, is among the highest redshift sources 
in our sample. It is much fainter in the H band lacking any 
morphological information. 


Figure 2 shows the spectral energy distributions of the 
two candidates shown in figure 1 together with their galaxy 
spectral fits and the derived photometric redshifts. It is to 
note how faint are the selected AGNs in the optical bands. 
The SEDs of all the selected AGN candidates are shown 
in the Appendix. This outlines the importance of the joint 
H band, X-ray detect ions to find out fai nt AGNs in mul¬ 
tiwavelength surveys (Fiore et al. (2012)). Figure 3 shows 
the X-ray over H-band hux ratio as a function of the H band 
magnitude for our GANDELS GOODS-S sample of AGN 
candidates. The straight line in the figure represents the 
locus of constant X-ray flux corresponding to the average 
detection threshold. Looking at the X/H ratio distribution 
derived at different H magnitudes it appears that only faint 
NIR z > 4 candidates which are relatively bright in the X- 
ray band can be detected in the Ghandra 4Ms field. In this 
respect our sample is biased against relatively X-ray faint 


AGNs, a feature that has to be taken into account in esti¬ 
mating the faint end of the AGN UV luminosity function 
as computed in the next section. We note again that our 
sample is coming from a well defined magnitude selected 
galaxy sample in the UV rest-frame band. The X-ray de¬ 
tection is only used in this context to reveal the presence 
of AGN activity. 


4. AGN UV luminosity functions 

The estimate of the faint end of the AGN luminos¬ 
ity function in the UV rest frame has been performed 
with an extended versio n of the standard 1/Umax al¬ 
gorithm (Schmidt (1968)1 where different regions in the 
GOODS-S hefd with different magnitude limits were com- 
bine d together in the volu me estimate of each object 
(e.g. |Avni & Bahcall (1980)|). Indeed, for each object and 
for each j-th region shown in table 1, an effective vol¬ 
ume Vmax{j) is computed. For a given redshift interval 
{ziow,Zup), these volumes are confined between ziow and 
ziim{j), the latter being defined as the minimum between 
Zup and the maximum redshift at which the object could 
have been observed within the magnitude limit of the j — th 
region. The galaxy number density z) in a given 

(Az, AM) bin can be computed as follows: 


(j){M,z) = 


1 


AM 


E 

2=1 






E 

dz 


dz 


( 1 ) 


where u}{j) is the area in units of steradians corresponding 
to the region j shown in table 1 , n is the number of objects 
in the bin and dV /dz is the comoving vol ume element per 
steradian (see e.g. Salimbeni et al. (2008) I. 

Rest frame UV 1450 A absolute magnitudes M 1450 were 
computed from apparent magnitudes in the filters closer to 
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A(A) A(A) 


Fig. 2. Spectral energy distributions of the objects as in figure 1 



Fig. 3. X/H flux ratio as a function of H magnitude for the GANDELS GOODS-S sample. The X-ray flux limit adopted 
in the Ghandra image is also shown together with the H magnitude limit H = 27. 


this rest frame wavelength for each redshift. For example, 
for AGNs at z ~ 4.5 the UV rest-frame absolute magni¬ 
tudes have been directly derived from the observed HST 
1775 magnitudes. The absolute magnitudes M 1450 used for 
the evaluation of the luminosity function are shown for each 
object in Table 2. 

The resulting volume densities as a function of M 1450 , z 
have then been corrected for various sources of incomplete¬ 
ness. The first correction is due to the gradual H band bend¬ 
ing of the galaxy counts which corresponds e.g. to a 50% 
drop at = 27 in area #4 of the present sample. In this 
respect, a correction factor has been included as a function 
of area and H magnitude in the volume estimate of each ob¬ 


ject. The second correction is due to the fact that, for each 
H magnitude, only relatively bright X-ray AGN candidates 
(with fx > 1.5 X 10“^^ erg s“^ cm“^) can be detected in 
our z > 4 parent sample (see Figure 3). This corresponds 
to the loss of AGN candidates with e.g. Fx/iX/n) ^ 0.1 
at H ^ 27. The incompleteness fraction is derived from 
the same X/Fl distribution observed above the X-ray flux 
threshold. The objects have been weighted considering that 
their H distribution in Figure 3 is affected by the faint slope 
of the UV rest-frame luminosity function and by the incom¬ 
pleteness in the H band counts. We adopted a faint end 
slope -1.5 and count incompleteness as derived in section 2. 
The two effects tend to compensate and the final correction 
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-19 -20 -21 -22 -23 -24 -25 -26 -27 -28 


Fig. 4. UV 1450 A AGN luminosity functions in various redshift bins. Different symbols represent different surveys as 
explained in the figure box. Open squares in the highest redshift interval (bottom panel) represent LF bins derived from 
AGN candidates with more uncertain photometric redshifts (see Section 6.1) 


is not sensitive to slope changes of ^ 0.2. Finally, correction 
factors to volume densities typically of the order of 10 - 20 % 
have been applied to volume densities after taking into ac¬ 
count spatial fluctuations in the X-ray flux limits for each 
position of the X-ray detection. 


Given the very poor statistics we only evaluated the 
Poisson contr ibution to the er rors in each LF bin adopting 
the recipe by Gehrels (1986) valid also for small numbers. 
In this respect the errors shown in the figure represent lower 
limits respect to the true values. Indeed for example cosmic 
variance can play a significant role in the present deep pen¬ 
cil beam survey increasing the uncertainties in the derived 
volume densities. 


The resulting luminosity functions are shown in figure 
4 and table 3 for M 1450 < —18.5 where the width and posi¬ 
tion of the bins have been selected to allow where possible 
a homogeneous statistics among the bins. It is clear that 
our faint AGN candidates are sampling the faint end of the 
AGN UV luminosity function in the absolute magnitude 
range —18.5 > M 1450 > —22.5 and with densities at the 


faint end between 10“^ ^ ^ 10“'* Mpc“^ mag“*. Thus 

the presence of even few faint high z AGN candidates in 
the GOODS-S field changes appreciably the abundance of 
the AGN population at z = 4 — 6.5 and its cosmological 
significance especially in the context of the cosmic reion¬ 
ization. 

In fact to provide a first estimate of the total AGN 
emissivity at z > 4 a shape of the luminosity function 
has been derived in each redshift bin connecting the vol¬ 
ume densities estimated from our sample with that of the 
brightest high z SLOAN QSO sample where selection ef¬ 
fects respect to the morphological appearance and X-ray 
properties are thought to be small. In other words, no X-ray 
QSOs with strong absorption in the rest-frame optical/UV 
are expected at the brightest magnitudes where the Sloan 
sample should be representative of the overall AGN pop¬ 
ulation. We adopted a double power-law shape to connect 
the two samples in each redshift bin of the type 


]^Q0.4(Mb„alc-M)(/3-l) _|_ X00-4(*fi.reofc-M)(7-l) 


( 2 ) 
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Table 3. AGN luminosity functions, emissivities and photoionization rates^ 


A2 

Mi450 

log (pobs 

log (pcorr 

No (,3 

// 

7 

"^^break 

log 0’ 

£1450 

£912 

r 

4 - 4.5 





1.52 

3.13 

-23.2 

-5.2 

11.45 

6.8 

0.62 


-19 

- 4.9 

-4.3 

1 









-20 

- 4.7 

-4.6 

3 









-21 

- 4.9 

-4.9 

2 








4.5 - 5 





1.81 

3.14 

-23.6 

-5.7 

9.90 

5.92 

0.43 


-19 

- 4.5 

-4.1 

3 









-20 

- 4.9 

-4.6 

2 









-21 

- 4.9 

-4.8 

2 









-22.5 

-5.3 

-5.3 

1 








5-6.5 





1.66^ 

3.35 

-23.4^ 

-5.8 

4.2 

2.5 

0.12 


-19 

-5.2 

-4.7 

1 









-20 

-5.1 

-4.7 

3 









-21 

-5.7 

-5.7 

1 









in units Mpc e in units of 10^“^ ergs s ^ Hz ^ Mpc ^ , F in units of 10 s (pcorr is the volume corrected for 
incompleteness in the X/H distribution. ^Value has been fixed to average in the two lower 2 intervals. 


The best fit parameters are shown in table 3. Note 
that in the highest redshift bin the faint end slope and 
the break magnitude have been fixed to the average val¬ 
ues in the two lower redshift intervals because of the lower 
statistics. In general however, it is clear from figure 4 that 
the break magnitude and the related volume density are 
not well constrained from the present composite sample. 
For this reason, the fitted volume density at the break is 
likely a lower limit and larger AGN surveys of intermedi¬ 
ate depth are needed to fix the point. In the lowest red¬ 
shift bin 4 < 2 < 4.5 we have also shown for compari¬ 
son the points of the luminosity function derived from a 
faint sample of type 1 QSOs with stellar appearance by 
Glikman et al. (2011) It is to note how the X-ray detection 
criterion used in the present paper is able to reach higher 
volume densities including AGNs with different morpho¬ 
logical and spectral properties. In the highest redshift bin 
5 < 2 < 6.5 we have also sho wn in Figure 4 the recent point 
by Kashikawa et al. (2014)| obtained from an optical spec- 

•4--v* 1 -V^T ry^*¥ r ttt l-i A ■ ' ^ y-J i y-1 r\ r\*r»y^ y-i y^ I y^ yn^-y*» 4-'v*y^'r-v-» 


5. AGN Hydrogen Ionizing Emissivity and 
Photoionization rate 

We now proceed to compute the contribution of our pre¬ 
dicted luminosity functions to the AGN hydrogen ionizing 


emissivity. This is given by 
iioniz) = {f)egi2 = 

if) J 4'{I^li50,z)Lii50 


0.44 


V 1200 / 


1.57 


dL 


1450 


( 3 ) 


where (/) is the average escape fraction of ionizing radi¬ 
ation from the AGN host galaxies and £912 is the emis¬ 
sivity produced by AGN activity. The shape of the AGN 
SED from A = 1450 A to A = 912 A has been rep¬ 
resented by a double power l a w with slopes ado pted 


after Schirber fc Bullock (2003) Telfer et al. (2002) and 


Vanden Berk et al. (2001) In the integral we have consid- 


troscopic survey where AGN candidates are selected from 
standard color selection, assuming a stellar appearance. We 
have added a vertical arrow to note that this value could 
represent a lower limit respect to combined NIR/X-ray sur¬ 
veys like that presented here. A decrease in volume density 
by a factor < 2 from 2 ^ 4.25 to 2 ^ 5.75 is present in 
our sample, smaller than found at higher luminositie s in 
UV and X-ray selected surveys (e.g. Vito et al. (2014) I. In 
summary, our new data give a first estimate of the overall 
shape of the UV AGN luminosity function in the absolute 
magnitude interval —18.5 > /U 1450 > —28. A flattening 
of the LF is needed to match the volume density of the 
brightest QSOs derived from the Sloan survey to that of 
our faint AGN population. The global behaviour of the lu¬ 
minosity function is used in the next section to derive the 
high redshift evolution of the UV emissivity. 


ered all the AGNs a factor 100 brighter or fainter than 
the break luminosity (corresponding to a faint limit of 
/Ui 45 o ~ —18). Gonsidering sources even 100 times fainter 
would increase the emissivity only by a factor ^ 10 % given 
the rather flat slope of the luminosity function. The uncer¬ 
tainties in the integrated emissivities are computed combin¬ 
ing the uncertainties in the normalization and bright end 
slope for the highest redshift bin since the faint end slope 
and break have been fixed. For the lower redshift bins the 
uncertainties are due to uncertainties on the normalization 
and faint end slope only (i.e. fixing the break luminosity 
and bright end slope). In this respect the widths of the er¬ 
ror bars in emissivity and photoionization rates should be 
considered as lower limits. It is to note that in the sam¬ 
pled redshift intervals the main contribution to the emis¬ 
sivity comes from AGNs near the predicted break magni¬ 
tude Mi 45 o —23.5. It is just at this magnitude where 
volume densities shoub be estimated with greater accuracy 
but this requires surveys of comparable depth extended on 
larger areas. 

In figure 5 and Table 3 the Lyman continuum ionizing 
emissivities £24 in units of 10 ^"^ erg s“^ Hz“^ Mpc“^ derived 
from the estimated AGN luminosity functions are shown as 
a function of redshift. We assume (/) = 1 as a reference 
value for the ionizing escape fraction as observed in most 
bright AGNs. 

The prediction of the model proposed by 
Giallongo et al. (2012)| is also 


is also shown in figure 5 for 
comparison. The model is based on the standard ACDM 
scenario for galaxy formation and evolution already 


described in Menci et al. (2005) and includes a physical 
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description of the Super Massive Black Hole growth 
at the center of galaxies as mainly due to the in- 
flow of cold gas destabilized by minor merger events 


(Cavaliere & Vittorini (2000) Menci et al. (2008) I. The 


model also includes a self consistent description of the 
ionizing emissivity of AGNs based on the blast-wave 
feedback mechanism ([Cavalie re, Lapi, fc Menci (2002)[ 


Lapi, Cavaliere, &: Menci (2005) f The blast-wave pro- 

duced in the interstellar medium sweeps out the gas 
allowing ionizing UV photons free to escape from the AGN 
host galaxy into the IGM. The model prediction tends 
to overestimate the observed emissivities by factors two 
to five as the redshift increases from 2:~4 toz^6 .5. 
In the same figure 5 the AGN model prediction by 


Haardt fc Madau (2012)| is also shown for comparison. 
The latter model underestimates by a larger amount the 
observed values up to a factor 20 at 2 ; ~ 6. Overall the 
new data are in between the two models, not far from the 


predictions by the Giallongo et al. (2012) model, and show 
a gradual decline of the ionized emissivity in the probed 
redshift interval 2 ; = 4 — 6.5. To look whether this decline 
is consistent with the evolution of the IGM ionization level 
we computed the hydrogen photoionization rate predicted 
from our emissivity and compared it with that inferred 
from the analysis of the Lyman-a forest in high redshift 
QSO spectra. 

The photoionization rate per hydrogen atom 

r_i 2 in units of 10“^^ s“^ was computed as in 


Giallongo et al. (2012) including the recent update on 
the mean free path (nifp) of ioni zing photons in the IGM 
derived by Worseck et al. (2014) 


r ^ ^ nr *^24(^:) f Al ^ f I + Z 

^ ' 3-1- I auv I v65Mpc/ \ 4.5 


3-r) 


where auv = —1.57 (Schirber & Bullock (2003) I. The 
value ai, = —1.57 adopted for 475 < A < 912 A is consistent 
with the recent evaluation with the COS spectrograph at 


HST by Shull et al. (2012) and Stevans et al. (2014)| a^ = 
— 1.41. The power-law index — ry describes the decrease in 
redshift of the proper mean free path (mfp) Al of ionizing 
photons in the IGM due to the increase in redshift of the 
Lyman Limit absorption systems. We adopt 77 = 5.4 and 
th e normalization of the m fp to 65 Mpc at z = 3.5 as found 
by Worseck et al. (2014) in the redshift range z = 4.4 —5.5. 
It is to note that the equation assumes that only ionizing 
sources within one absorption length contribute to the cos¬ 
mic reionization. The values derived from our GANDELS 
GOODS-S sample of AGN candidates show in figure 6 a 
gradual decline from z = 4toz = 6.5 although slightly 
steeper than that derived for the emissivity due to the in¬ 
crease of the average IGM absorption with increasing red¬ 
shift. 

Different data points derived in a mo del-dependent way 
from the IGM statistics observed in the Lyman-a forest 
of high z QSO spectra are also shown for comparison as 


Giallongo et al. (2012) Some values depend on the UV 


background adopted in numerical simulations to repro- 


forest of QSO spectra ( 

Faucher-Giguere et al. (2008) 

Bolton et al. (2005) 

Becker & Bolton (2013) 

b'a.n et al. (2006) Wyithe & Bolton (2011)). Other values 


the highest redshifts z ^ 5 — 6 by Calverley et al. (2011) 


Different methods produce the significant scatter shown 
in figure 6 despite the data are scaled to the same IGM 


temperature-density relation as in Giallongo et al. (2012) 
This is due to uncertainties on modelling the gas den- 
sity and temperature distributions in the IGM at very 
low densities ( Miralda-Escude, Haehnelt, fc Rees (2000)[ ) 
and within the proxim ity region surrounding each QSO 
(Calverley et al. (2011)). 

It is important to note from figure 6 that the photoion¬ 
ization rates at z > 4 derived by our AGN luminosity 
functions are consistent with the photoionization levels of 
the IGM derived from the Lyman-a statistics in the same 
redshift range. This would imply a reionization occurring 
at redshifts not much larger than z ~ 7 since moderately 
bright ionizing sources with a maximum comoving density 
of ^ 2 X 10“® Mpc“^ at z ~ 6 have a typical separa¬ 
tion of ^ 3 proper Mpc, only a factor 2 smaller respect 
to the derived prop er mean free path of io nizing photons 
in the IGM at z ^ 6 Worseck et al. (2014) The latter sce¬ 
nario would also be consistent with the large decline of 
strong Lyman-a emitting galaxies observed between z ^ 6 
and z ^ 7 which would suggest a large neutral hydrogen 
fraction in the IGM and a patchy reionization process al 
ready at z 7 dPentericci et al. (2011)||Treu et al. (2013) 


Pentericci et al. (2014) Tilvi et al. (2di4]|T 


6. Discussion 

In the context of this preliminary analysis we should dis¬ 
cuss some caveats which could significantly reduce our ex¬ 
pectations about the contribution to the reionization by 
high redshift AGNs of intermediate-low luminosities. First 
of all some redshift estimates for our candidates could be 
wrong, despite the average high level of accuracy, reduc¬ 
ing the number of high z AGNs. The second issue concerns 
the assumption about the predominance of the AGN UV 
spectral energy distribution even in cases where the optical 
spectra appear affected by significant stellar contribution 
or look very steep. A further issue concerns the assumption 
of a high escape fraction of ionizing photons for the global 
AGN population, including the X-ray absorbed AGN frac¬ 
tion. Finally we discuss the implication that an ionizing 
AGN population at very high redshift would have for an 
early Hell reionization of the IGM. 


6.1. Redshift reliability 

Goncerning the reliability of the derived redshift distribu¬ 
tion in our sample we should note that the evaluation of 
the photometric redshifts becomes progressively more un¬ 
certain for fainter sources with featureless SEDs. For this 
reason we have shown in the Appendix the probability red¬ 
shift distributions PDF(z) for our candidates. In this con¬ 
text, as already stated, the estimate of photometric red¬ 
shifts for sources at z > 4 mainly relies on the statistical 
significance of the Lyman a forest (< 1216 A) and Lyman 
break at 912 A rest frame wavelength. In fact, the expected 
escaping Lyman continuum emission from the sources, even 
assuming (/) ~ 1, would be strongly depressed by IGM 
absorption causing a flux dropout. As a consequence, the 
redshift uncertainty at z = 5 — 6 is related to the flux 
dropouts near the Lyman-a and Lyman edge almost inde¬ 
pendently of the assumed galaxy or AGN spectral library. 
Our assumption is corroborated by the good agreement we 
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found between the 5 available spectroscopic redshifts and 
the photometric estimates of the GANDELS catalog. It is 
clear however that in cases where the spectrum is partic¬ 
ularly steep the evidence for the presence of any Lyman 
break weakens, increasing consequently the uncertainty in 
the redshift estimate. This is particularly true for some of 
the z > 5 objects. Indeed we have already excluded from 
the LF analysis the object 29323 because of its too peculiar 
SED and PDF(z). Thus five sources have been used for the 
LF estimate in the highest redshift bin, among them 20765, 
28476 and 33160 have the most uncertain redshift estimate 
as shown in Figure A.l. To get a rough estimate of the un¬ 
certainties involved in the derived volume densities of faint 
AGNs at z > 5 we have repeated the estimate of the lumi¬ 
nosity function at the highest redshift bin excluding these 
sources. Two of them are the only sources in the faintest 
and brightest LF bins (M 1450 = —19 and M 1450 = —21) 
that would be removed. We have indicated with different 
symbols these uncertain LF bins. Two sources remain in 
the LF bin at M 1450 = — 20 slightly decreasing the aver¬ 
age volume densities by —0.1 to \og(j) = —4.8. This value 
would still be consistent with the two-power law extrapo¬ 
lation adopted to estimate the UV emissivity. Finally we 
note that some of these sources are relatively bright at 8 
micron but this by no means represents a prior against a 
high redshift solution for the estimated redshift. Indeed one 
of the best st udied X-ray absorbed f aint AGN in our cat¬ 
alog, 273 (e.g. 
observed by A^ 


Vanzella et al. ( 2008)|), w hich has also been 
MA (Gilli et ai. (2013)), has a robust spec¬ 
troscopic redshift at 2 = 4.76 and an 8 micron apparent AB 
magnitude of 21.5. Also the source 14800 at spectroscopic 
redshift z = 4.82 shows a relatively bright IRAG continuum 

at levels of 22.5. _ 

Very recently Weigel et al. (2015) have searched for z > 
5 AGNs in the same GUODS-S held using almost the same 
GANDELS dataset finding no convincing AGN candidates. 
The different result depends on their adopted procedure. 
They have looked for z > 5 sources starting from the 


Xue et al. (2011) X-ray selected catalog, looking for plau- 
sible optical drop-outs and/or photometric redshifts in the 
GANDELS images. Our sample is based on NIR selection in 
the H band and reaches fainter X-ray fluxes respect to the 


Xue et al. (2011) catalog. Moreover photometric redshifts 
in our sample have been obtained from the GANDELS mul¬ 
tiwavelength catalog derived from the B to the 8 micron 
SPITZER band that included careful SPITZER deblended 
photometry. As an example one source in our sample, (9713 
not included in the Xue et al. (2011) ) has a spectroscopic 
confirmation at z = 5.7. In our sample of 22 AGN can- 


didates only 2 are in common with the Xue et al. (2011) 
catalog with an estimated redshift z > 5. One of the two 
29323 (X156), has already been excluded from our analysis 
as stated above. The second object, 33160 (X85), appears 
at z ~ 6 although with a larger uncertainty and we have 
already checked in the previous paragraph the effects of its 
removal from the sample. 


6.2. UV spectral shape 

With few broad band filters often showing only upper lim¬ 
its it is impossible to disentangle the relative AGN and 
stellar contributions to the SED of our faint H = 24 — 27 
mag sources. This is why we do not include in our sample 
AGN candidates at z < 4. Indeed only at very high redshifts 


(z = 4—6) the photometric redshifts are mainly constrained 
by strong Lyman absorptions by the intergalactic medium 
intervening along the line of sight, almost independently of 
the intrinsic spectral shape assumed (AGN or galaxy tem¬ 
plate). If the sources are at very high redshifts, the observed 
optical bands (< 1 micron) give access to the rest frame UV 


SED shortward of 3000 A. The presence of X-ray detection 
in our sample at levels log Lx ^ 43 (erg s“^) implies prob¬ 
able AGN emission. As a consequence, somewhere in the 
UV part of the spectrum the AGN emission should emerge 
above the stellar spectrum. We assume that this happens 
at wavelengths just shortward of 1500 A. In particular we 
assume that our 1450 A luminosity is mainly due to AGN 
emission. Of course to verify this hypothesis a detailed spec¬ 
troscopic information would be required which is beyond 
the capability of the present instrumentation due to the 
faintness of our candidates. Nevertheless this assumption is 
supported by the distribution of the rest frame XJUV flux 
ratio in our sample (which at z^4 corresponds to the X/H 
ratio distribution shown in Figure 3). Our AGN candidates 
are confined in the X/H interval between about 0.01 and 
0.3-0.5 (in the same units of Figure 3) usually occupied 
by the AGN population in brighter AGN surveys as shown 
e.g. in Fiore et al. (201^ fo r the GOSMOS high redshift 
AGNs (|Givano et ai. (2011)). Also the few SLOAN QSOs 
at z = 6 show X/H fluxes which occupy the same region 
up to 0.2 — 0.3 although most of them have ratios below 0.1 
because they are very bright optically selected sources (and 
so with a significant X-ray weak fraction). If the rest frame 
UV fluxes were dominated by stellar emission this would 
imply an intrinsic AGN X/H flux ratios much (e.g. ^10 
times) greater than generally observed in brighter high red¬ 
shift AGNs. Of course both rest-frame X-ray and UV fluxes 
could be produced by stellar populations but this seems un¬ 
like for X-ray luminosities log Lx ^ 43 (erg s“^). By the 
way the suggestion provided by the SED models about the 
possible presence of Balmer breaks appears not compelling 
and not necessarily caused by stars. Rest-frame optical neb¬ 
ular emission lines (at z = 5 — 6 ) could mimick Balmer 
breaks whose contribution to broad band filters increases 


with (1 -I- z) as noted by Shaerer & de Barros (2009) 


Finally to check if the average spectral shape assumed to 
estimate the ionizing emissivity from the UV AGN luminos¬ 
ity function is appropriate for our sample, we have derived 
a first rough guess of the average power-law slope in the UV 
rest-frame interval 1300-1700 A, just around 1500 A. The 
average slopes have been derived from e.g. the V,I,z bands 
at z ^ 4 and Y,J bands at z ^ 6 . The average/median 
slope ay was found -1.4/-1.2 (excluding object 9323) al¬ 
though the dispersion is significant. In our paper we have 
adopted a flatter slope from 1450 down to 1200 A (-0.44) 
and a steeper one shortward of 1200 A (-1.57). Assuming 
the derived average slope -1.4 (instead of -0.44) just short¬ 
ward of 1450 A reduces the photoionization rates by about 
10 percent only. To halve the derived photoionization rates 
in Fig .6 we should adopt a much steeper -2.5 slope for all the 
AGNs in the luminosity function. The average UV slope de¬ 
rived in our sample,-!.4, appears thus consistent with that 
(-1.57) adopted in the paper for the overall AGN population 
from i\Li 45 o = —18 to Mi 45 o = —28. 


In summary, although we can not reliably quantify the 
stellar contribution to the optical/UV emission we assume 
that in our X-ray sources the SED shortward of the rest- 
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frame 1450 A is dominated by AGN emission. This assump¬ 
tion appears consistent with the typical AGN X/H (X/UV 
rest-frame) flux ratio found in our sample, which would im¬ 
ply an average UV-X spectral slope similar to that found 
in AGNs. 


6.3. Ionizing escape fraction 

One of the main uncertainties concerns the assumption 
about the escape fraction of UV photons from AGNs. In 
fact, we assume that all ionizing photons are able to es¬ 
cape and ionize the surrounding IGM. This is usually al¬ 
most true for bright type 1 QSOs where the optical depth 
is Tii < 1 at the Lyman limit position corresponding to 
the systemic redshift of the QSOs. Any appreciable absorp¬ 
tion is present at shorter wavelengt hs due to the interven - 
ing IGM along the line-of-sight (e.g. |Sargent et al. (1989) I. 
Our faint sample of AGN candidates lacks spectroscopic 
information except for the 5 cases shown in table 2 and 
could include a significant fraction of absorbed AGNs. Hints 
of excess absorption have been derived in the X-ray spec¬ 
tral analysis of high redshift X-ray selected AGNs but the 
reliability of any measure is hampered by the fact that 
the photoelectric cutoff of the aborption goes out of the 
spectral region sampled by the Chandra observations al¬ 


ready at z > 3 (see e.g. Vito et al. (2013) I. Moreover to 


(Hasinger (2008) Merloni et al. (2014) I although there is 


no tight correlation between X-ray absorption and the pres¬ 
ence of broad emission lines since some X-ray absorbed 
AGNs show broad emission lines in their rest-frame UV 


spectra, as outlined by Fiore et al. (2012) Moreover recent 
COS spectral analysis of FUV spectra of low z AGNs of 
different luminosities (log L[/y = 44 — 46.5) and spectral 
typ es (Sy 1, Sy 1.2, Sy 1 .5, QSOs) by Shull et al. (2012) 


and Stevans et al. (2014) show no continuum edges at 912 
A. In this complex situation what we are assuming is that 
the presence of any AGN activity in high z, rest-frame UV 
selected galaxies - activity certified by X-ray emission at 
levels log Lx > 43 - is able to provide the needed mechani¬ 
cal and radiative feedback mechanisms to allow ionizing UV 
photons free to escape outside the host galaxies into the sur¬ 
rounding IGM with average escape fractions in the sample 
not much smaller than 1 (e.g. in the range (/) ^ 0.7— 1). It 
is to note that a reduction of the expected photoionization 
rate by e.g. 30% respect to that shown in Figure 6 , due to an 
assumed average escape fraction (/) ^ 0.7, would be com¬ 
pensated by the re-emission of io nizing photons by the IGM 
due to radiative recombination (Haardt & Madau (1996)) 
that in a framework of reionization by the AGN population 
could contribute by a factor ^ 1.25 to the cosmological 
photoionization rate at z > 5 (Haardt, private communica¬ 
tion). 


allow a first estimate of the hydrogen column density with 
few photon counts, a fixed slope (typically F = 1 . 8 ) of 
the intrinsic X-ray power-law is often adopted and a solar 
metallicity is always assumed despite observational hints 
of supersolar metallicitie s in the black hole neighbourhood. 
In fact flatter F values (|Vito et al. (2013)|) and supersolar 
metallicities could both significantly reduce the estimated 
hydrogen column density. Some statistics on the relative 
fraction of obscured/unobscured QSOs of intermediate X- 
ray luminosity (43 < logL^ ^ 44) shows that the frac- 
tion of obscured AGNs at high redshift s could be ^ 50% 


Feedback at various intensity levels is indeed emerg¬ 
ing from the most recent observations as an ubiquitous ac- 
tivity in AGNs, inde pendently of their spectral type (e.g. 
Gicone et al. (2014) and reference therein). Due to the feed- 
back mechanisms in the AGN host galaxies some escaping 
ionizing photons could also come from young stars in the 
hosts. For a given AGN/stellar flux ratio at 1500 A, the stel¬ 
lar ionizing spectrum at A < 912 A is expected to be softer 
than shown by the AGN due to intrinsic Lyman continuum 
absorption in the stellar atmospheres. The latter absorption 
is essentially unknown and expected values depend on the 
spectral synthesis models adopted, on the age and metal¬ 
licity assumed for the stellar populations in star forming 
galaxies. Typical values for the F 1500 /A 912 flux ratio, which 
is a measure of the intrinsic stellar absorption, are 3 for 
normal star forming galaxies. However for our AGN host 
galaxies at z > 4 the age and metallicity of the associated 
starburst activity could be < 10 ^ yr and ^ O. 2 Z 0 , respec¬ 
tively, implying intrinsic F 1500 /F 912 ratios ^ 2 or smaller. 
These values could be not far from the ratio 1.7 derived 
from our adopted AGN spectral shape. In other words the 
ionizing UV spectra from a very young population in AGN 
host galaxies could be almost as hard as that of the AGNs 
at wavelengths just shortward of the Lyman edge at 912 A. 
Of course with decreasing wavelengths shortward of 912 A 
the AGN and stellar spectral shapes are progressively more 
and more different and at the Hell edge (228 A) only the 
AGN SEDs can contribute to the Hell reionization. 


6.4. Hell reionization 

For this reason, in our proposed scenario for the cosmic 
reionization, the hardness of the ionizing UV spectrum 
provided by the AGN population would produce an early 
significant reionization of the Hell content in the IGM. 
In fact, there is some consensus that Hell reionization 
could be completed later in the universe, at z ^ 3, re¬ 
spect to hydrogen, supporting a scenario where galaxies 
with much softer UV spectra dominate the ionization of 
hydrogen at very high redshifts until AGNs become so 
numerous at z 3 to provide the required hard pho- 
tons to reionize Hell. The issue has been already discussed 


Giallongo et al. (2012) in the framework of an AGN 


dominated U V background where the Helll volume fill¬ 
ing factor was computed as a function of decreasing red¬ 
shift under various assumptions. In particular it was rec¬ 
ognized that the presence of optically thick Hell clouds 
between AGNs reduce the mean free path of Hell ioniz- 


ing photons, according to Bolton, Oh, &: Furlanetto (2009) 
allowing a gradual increase of the Hell reionization as 
the redshift decreases up to a full reionization at z ^ 3. 
Of course a detailed evolutionary history of the hydro¬ 
gen and helium reionization requires 3D hydrodynami- 
cal models however the simplified model suggests that a 
modest degree of inhomogeneity in the Hell distribution 
could delay its reionization history. In this context, recent 
direct observations of Hell spectr al regions in 17 bright 
QSOs at z ^ 3 with HST/COS by Worseck et al. (2014b) 
seem to suggest a mild evolution of the Hell reioniza- 
tion with increasing redshift allowing the possibility to 
get the bulk of the Hell reionization in the universe at 
z > 4. On the other hand rec ent detailed hydrodynam- 
ical simulations of the IGM (Gompostella et al. (2014) 
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® Faucher-Giguere et al. 2008 ^ Kirkman et al. 2005 

BSolton et al. 2005 ^ Calverly et al. 2011 

□ Becker &c Bolton 2013 ^ Wyithe &c Bolton 2010 




Fig. 5. Cosmic ionizing emissivity by AGNs as a function of 
redshift assuming (/) = 1. Black squares are from our sam¬ 
ple. Red continuous curve is from the Haardt and Madau 
(2012) model. The long dashed curve is from the Giallongo 
et al. ( 2012 ) model. 


Puchwein et al 
of the hydrogen 


. (2014)1 ) 
n and H« 


enligthen the need for an increase 
hll photoionization rates respect to 


what predicted by the Haardt & Madau (2012) model at 
very high redshifts. These analyses if confirmed could sup¬ 
port a reionization scenario driven by AGNs. 


7. Summary 

Thanks to the completion of deep multiwavelength surveys 
it is now possible to apply specific selection procedures to 
reveal the very high redshift AGN population with intrinsic 
luminosities more typical of local Seyfert galaxies rather 
than bright Sloan QSOs. 

More specifically we have selected AGN candidates at 
z > 4 in the GANDELS GOODS-S field starting from the 
NIR H band selection of the galaxy parent sample, which 
corresponds to UV rest-frame selection, down to a very faint 
magnitude limit H < 27. We have used photometric and 
few spectroscopic redshifts derived from the GANDELS 
catalog and we have finally extracted from these high z 
galaxies the subsample of AGN candidates showing X-ray 
detection above Fx ^ 1-5 x 10“^^ erg cm“^ s“^ (0.5-2 
keV), correponding to a probability of spurious detection 
of 2 X 10“"^ in the Ghandra 4Msec GOODS-S field. 

We have made a preliminary estimate of the selection 
function aimed at deriving a correction for the sample in¬ 
completeness. We have derived average volume densities 
10 “® ^ ^ 10 “"* Mpc“^ mag“^ in the magnitude inter¬ 

val — 18.5 > Mi 45 o > —22.5. Thus the presence of even 


Fig. 6. Gosmic ionizing photoionization rate r_i 2 in units 
of 10“^^ s“^ produced by AGNs as a function of redshift as¬ 
suming (/) = 1. Black filled squares represent the predicted 
contribution by faint AGNs from the GOODS-S sample. 
Other small red symbols are the values inferred from the 
ionization status of the IGM as derived from the Lyman a 
forest analysis in high z QSO spectra. 


few faint high z AGN candidates in the GOODS-S field 
increases appreciably the cosmological significance of the 
AGN population at z = 4 — 6.5 especially in the context of 
the cosmic reionization. 

To make a first estimate of the ionizing UV emissivity 
of the AGN population at these high z we have derived 
for the first time the shape of the UV luminosity function 
from Mi 45 o ^ —18.5 to M 1450 ~ —28 adding the bright 
QSO samples from the Sloan survey which is less affected 
by significant incompleteness. The UV luminosity function 
is well represented by a double power-law although the po¬ 
sition of the break is highly uncertain. The values obtained 
from a simple fit are used to derive an ionizing emissivity 
of the order of 624 ^ 6.5 at z = 4 — 5 and 624 ^ 2.5 at 
z ~ 6 , assuming an average AGN UV spectral shape and 
an escape fraction of ionizing UV photons (/) = 1. 

The derived emissivity has then been used to estimate 
the hydrogen photoionization rate in the same redshift in¬ 
terval which is of the order ofr_i 2 ~ 0.1 at z ~ 6 , consis¬ 
tent with the value required to keep highly ionized the inter- 
galactic medium observed in the Lyman-a forest of bright 
QSO spectra at the same redshift. This scenario points to¬ 
ward a late reionization epoch at redshifts not much higher 
than z = 7 as suggested by the rapid decrease in the number 
of Lyman alpha emitters and by the recent lower thomp- 
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son OTtical depth suggested by recent preliminary Planck 
datcfl 

We have discussed some caveats which could reduce the 
volume density and/or ionizing emissivity of the high red- 
shift AGN population described in the present work. More 
specifically the large uncertainties involved in photometric 
redshifts for sources with steep, almost featureless, spec¬ 
tral energy distribution could spuriously enhance the to¬ 
tal number of z > 5 AGN candidates in the GOODS-S 
field. However, after careful evaluation of this bias, we have 
shown that possible changes in the AGN volume densities 
provide values still consistent with the adopted shape of the 
UV luminosity function in the redshift interval 5 < z < 6.5. 
A further important caveat is related to the escape frac¬ 
tion of ionizing photons from the AGN host galaxies which 
should be on average (/) > 0.5 to ensure enough ionizing 
photons in the universe. At present it is not clear if fainter 
AGNs ionize their neighborhood as well as the brighter 
QSO population does. 

We finally emphasize that, given the large uncertainty in 
the present data, the aim of this preliminary study is not fo¬ 
cused on deriving a reliable shape of the AGN UV luminos¬ 
ity function but is rather devoted to establish if faint AGNs 
can be definitely ruled out as main contributors to the cos¬ 
mic reionization, as generally thought. The main outcome 
of the present preliminary analysis goes in the opposite di¬ 
rection. Given the large uncertainties involved there is room 
for a significant contribution to reionization by the z > 4 
AGN population under reasonable physical assumptions. 
Looking for faint NIR selected high-z galaxies with even 
marginal X-ray detection in the deepest available NIR and 
X-ray images can be a promising technique to reveal a large 
number of Lyman Gontinuum emitting sources. Applying 
the same selection strategy and spectroscopic follow up to 
other fields with similar characteristics like e.g. GOODS-N 
and GOSMOS will provide a composite larger sample and 
consequently a more reliable estimate of the AGN UV lu¬ 
minosity function and photoionization rates up to redshift 
z > 6 where the cosmic hydrogen reionization is still in 
action. 
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Appendix A: Spectral energy distributions and 
redshift probability distribution functions 

In this section we show the spectral energy distributions 
(SEDs) and the probability distribution functions (PDFs) 
of the photometric red shift estimates for ea ch AGN can¬ 
didate derived from the Dahlen el al. (2013) analysis. Most 
of the candidates in figure Al show PDh's confined at z > 4 
with only small wings at z < 4. The most uncertain red- 
shifts being for the objects 20765, 28476, 29323, 33160 al¬ 
though the larger PDFs distributions appear mostly dis¬ 
tributed at z > 4. It is to note that the latter three ob¬ 
jects at z > 6 are among the most uncertain due to com¬ 
bination of power-law shape of their SEDs and faintness of 
the sources. Note also that for some of these objects, the 
bayesian photometric redshift does not correspond to the 
peak of the PDF but it is the result of a weighted average 
around the maximum. 


Appendix B: X-ray position accuracy 

In this section we show an example of the relative position 
accuracy between the X-ray and H band image. We have 
selected the object 28476 already shown in Figure 1 which 
appears located near two brighter H band sources. In Figure 
B.l we also show the position accuracy of a bright source 
on the left bottom corner which results comparable to that 
of 28476. In general the position accuracy between the H 
and X-ray sources is < 1 arcsec for all the AGN candidates 
of our sample. They have off-axis Ghandra positions < 9 
arcmin. 
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Fig. A.l. Spectral energy distributions and probability dis¬ 
tribution functions of photometric redshifts normalized to 
1 at the peak for all the AGN candidates shown in Table 
2. 1-(T upper limits are also shown as downward arrows. 



Fig. B.l. H band image around the AGN candidate 28476 
(upper right) shown in Figure 1. The size is 30 arcsec. Gyan 
and green contours are X-ray detections in the 0.8-4 keV 
and 0.5-2 keV bands, respectively. Position accuracies are 
within 1 arcsec. 













































